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A statistically-selected Chandra sample of 20 galaxy clusters - 1. 
Temperature and cooling time profiles 
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ABSTRACT 

We present an analysis of 20 galaxy clusters observed with the Chandra X-ray satellite, 
focussing on the temperature structure of the intracluster medium and the cooling time of the 
gas. Our sample is drawn from a flux-limited catalogue but excludes the Fornax, Coma and 
Centaurus clusters, owing to their large angular size compared to the Chandra field-of-view. 
We describe a quantitative measure of the impact of central cooling, and find that the sample 
comprises 9 clusters possessing cool cores and 1 1 without. The properties of these two types 
differ markedly, but there is a high degree of uniformity amongst the cool core clusters, which 
obey a nearly universal radial scaling in temperature of the form T oc r~° ' 4 , within the core. 
This uniformity persists in the gas cooling time, which varies more strongly with radius in 
cool core clusters (t coo i oc r~ 1-3 ), reaching t coo \ < lGyr in all cases, although surprisingly 
low central cooling times (< 5Gyr) are found in many of the non-cool core systems. The 
scatter between the cooling time profiles of all the clusters is found to be remarkably small, 
implying a universal form for the cooling time of gas at a given physical radius in virialized 
systems, in agreement with recent previous work. Our results favour cluster merging as the 
primary factor in preventing the formation of cool cores. 

Key words: galaxies: clusters: general - intergalactic medium - X-rays: galaxies clusters 



1 INTRODUCTION 

Clusters of galaxies are uniquely suitable environments for study- 
ing galaxy formation and evolution, as well as excellent cosmologi- 
cal tools for measuring the fundamental parameters of the Universe. 
A study of the scaling properties of clusters, which span roughly 
two decades in mass, is one of the best means of tackling both these 
objectives, and many analyses have been conducted on this basis, 
establishing broad trends and quantifying departures from simple 
self-similarity (e.g. Edge & Stewart 1991b; Markevitch et al. 1998; 
iHelsdo n & Pon maj200d:lFukazawa. Makishima & Ohashil2004l) . 

Increasingly, however, attention is shifting towards a more 
complete understanding of cluster physics, focussing in particular 
on the scatter in these scaling relations. With current generation 
telescopes, we are reaching the point where the precision of scal- 
ing relations is limited by systematic effects rather than measure- 
ment errors. Further progress requires a better knowledge of clus- 
ter physics to understand the sources of intrinsic scatter in these 
relations. Not only does this hold great potential for understanding 
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galaxy evolution, it also promises to hone cluster observables as 
proxies for mass, which will be vital in extracting precise cosmo- 
logical constraints from forthcoming large cluster surveys. 

Since the cooling time of gas in the centres of most clusters 
is short compared to their age, radiative cooling would be expected 
to play an important role in cluster physics. However, many clus- 
ters do not show strong signs of cooling, possibly as a result of 
disruption caused by merging and associated shock heating of the 
gas, or else non-gravitational energy input (e.g. from AGN), or even 
heat redistribution via conduction or mixing. Moreover, even where 
cooling dominates, its progress is apparently impeded (see the re- 
cent review by P eterson & F abian 2006, and references therein), by 
processes which may well relate to these same phenomena. What- 
ever the explanation, it is clear that clusters can be separated into 
two distinct categories, accordin g to the presenc e or absence of a 
'cool core' iPeres et aljfl998l iBauer et al .0 2005 ). where radiative 
heat losses have significantly lowered the central temperature of the 
gas. It is the combination of cool core and non cool core systems 
in the cluster population which accounts for much of the intrinsic 
scatter in scaling relations. 

While cool core clusters have tended to be relatively well- 
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studied, non-cool core clusters have received less attention in scal- 
ing relation analyses. Consequently there is great potential for ex- 
ploring cluster physics by comparing the properties of the two types 
of clusters. This aim of this work is to investigate precisely this 
comparison, by focussing on a systematically selected sample of 
clusters, using data from the Chandra satellite, to allow a high- 
resolution study of the core properties of the intracluster medium 
(ICM). 

Previous detailed studies of the temperature structure of the 
hot gas in clusters have generally concentra ted on the m ost relaxed 
(i.e. morphologically regular) clusters (e.g . ISanderson et alj|2003t 
Piffar etti et alj|200a IVikhlinin et alj|2005l) . which frequently pos- 
sess cool cores. In this work, we have compiled a statistically se- 
lected sample of 20 clusters, to provide a more rigorous selection 
process. This naturally includes a larger number of non-cool core 
and perhaps less-relaxed clusters, which are more representative 
of the cluster population as a whole. We have subjected these 20 
systems to a non-parametric analysis, so as to avoid any model- 
dependent biases in the comparison of cool core (CC) and non-CC 
systems. Furthermore, in restricting the analysis to X-ray data taken 
with the Chandra satellite, we are able to exploit the unrivalled spa- 
tial resolution of this telescope, to permit a detailed examination of 
the inner regions of the ICM, where the impact of radiative cooling 
is most pronounced. 

Throughout this paper we adopt the following cosmological 
parameters; Ho = 70 km s" 1 Mpc~ , fl m = 0.3 and C2a = 0.7. 
Throughout our spectral analysis we h ave used XSPEC 11.3.1, in - 
corporating the solar abundance table o flGrevesse & SauvalUl998l) . 
which is different from the default abundance table. Typically this 
results in larger Fe abundances, by a factor of ~1.4. All errors are 
la, unless otherwise stated. 



2 SAMPLE SELECTION AND PROPERTIES 

The sample comprises the 20 highe st flux clust e rs draw n from the 
63 cluster, flux-limited sample of Ikebe et al. 1 2002), excluding 
those objects with extremely large angular sizes (the Coma, For- 
nax and Centaurus clusters), which are difficult to observe with 
Chandra owing to its limited field of view. So me basic p roper- 
ties of this sample are listed in Table Q The Ikeb e et alj flux- 
limited sample was c onstructed from the HIFLUGS sample of 
iReiprich & Bohringed 120021) . additionally selecting clusters lying 
above an absolute galactic latitude of 20 degrees and located out- 
side of the Magellanic Clouds and the Virgo Cluster regions. 

Fig. Q shows the cluster redshift as a function of mean X-ray 
temperature, excluding any cool core (see Section l4~TY The temper- 
ature errors are generally too small to be seen, and have been omit- 
ted for clarity. The point style differentiates between those clusters 
possessing a significant cool core and those without. The cool core 
classification scheme is described in Section l4~2l 



3 DATA REDUCTION 

The data analysis and reduction were performed with version 3.2.2 
of the standard software — Chandra Interactive Analysis of Ob- 
servations (CIAO 1 ), incorporating CALDB version 3.1. For all the 
observations a new level 2 events file was generated from the level 
1 events file downloaded from the Chandra archive. This ensures 

1 http://cxc.harvard.edu/ciao/ 
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Figure 1. Cluster mean temperature within 0.1-0.2 r$oo as a function of 
redshift, labelled according to the presence of a cool core. Temperature er- 
rors have been omitted for clarity. Redshifts are taken from the NASA Ex- 
tragalactic Database (NED). 

that the latest calibration information is applied uniformly to the 
data, irrespective of when the observations were made. 

The following procedure was followed for each observation 
dataset. According to which CCD chips were used in the obser- 
vation, a different light curve was extracted for CCD chips 5 and 
7 separately and the remaining, front-illuminated chips combined. 
The recommended criteria for energy extraction and time binning 
were used 2 . Flares were identified and excluded with the sigma 
clipping algorithm implemented in the CIAO task 'lc_clean', using 
the median light curve value to provide a robust initial estimate of 
the quiescent mean level. 

Cosmic ray events were identified and excluded using the 
CIAO task 'acis_run_hotpix', and those found in this manner were 
extracted separately and examined to check that photons from the 
cluster itself were not misidentified. A new level 1 events file 
was produced by reprocessing the resulting events file, to apply 
the latest gain file. Corrections were applied for the effects of 
charge transfer inefficiency (CTI) and time-dependent gain varia- 
tion, where necessary. Bad columns and hot pixels were then ex- 
cluded, and only events with Advanced Satellite for Cosmology 
and Astrophysics (ASCA) grades 0, 2, 3, 4, and 6 were retained. 
Subsequently, a new level 2 events file was generated by reprocess- 
ing this modified level 1 events data set. For those observations 
telemetered in very faint (VF) mode (indicated in column 3 of Ta- 
ble 0, the extra background event flagging and removal that this 
enables was performed in both the main and corresponding blank 
sky datasets. 

Owing to the close proximity of the clusters in this sample, 
emission from the target fills the entire Chandra field of view in 
most cases. It was therefore necessary to employ separate back- 
ground events files for each observation, using the Markevitch 

2 http://cxc.harvard.edu/contrib/maxim/bg 
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Name 


Obsid" 


Detector 6 


RA 
(J zuuu ) 


Dec. 
/nonn~i 


Redshift 


HI Column c 
t x iu cm ) 


Mean kT d 
^Ke v ) 


rsoo 


Notes' 2 


NGC 5044 


798 


S 


198.850 


-16.385 


0.008 


4.9 


i ig +003 
-0.02 


436 


- 


Abell 262 


2215 


S (VF) 


28.194 


36.152 


0.016 


5.4 


, 34 +0.06 
-0.06 


668 


- 


Abell 4038 


4188 


I(VF) 


356.928 


-28.143 


0.030 


1.6 


3 02 + - 17 
J - uz -0.17 


784 


F,T 


Abell 1060 


2220 


I(VF) 


159.181 


-27.526 


0.012 


4.8 


, ,0+0.07 
3 - ly -0.07 


812 


F,T 


Abell 1367 


514 


S 


176.191 


19.699 


0.022 


2.3 


3 n o +0.11 
i Z>i -0.10 


826 


T 


2A0335+096 


919 


S 


54.669 


9.970 


0.035 


17.9 


3 47+0.14 
3 - H ' -0.12 


856 


- 


Abell 2147 


3211 


I(VF) 


240.567 


15.963 


0.035 


3.4 


4 17 +o.i4 

*- LI -0.14 


961 


F,T 


Abell 2199 


497 


S 


247.160 


39.551 


0.030 


0.9 


4 36 +0 - 13 
^■ 30 -0.11 


989 


- 


Abell 496 


931 


s 


68.408 


-13.262 


0.033 


4.6 


5.M ±8:5! 


1200 


- 


Abell 2256 


1386 


I 


256.041 


78.648 


0.058 


4.1 


6.08± -3 3 ? 


1220 


T 


Abell 1795 


493 


S (VF) 


207.219 


26.590 


0.062 


1.2 


6 30 + 17 

OJU -0.17 


1250 


- 


Abell 3558 


1646 


S (VF) 


201.987 


-31.496 


0.048 


3.9 


6 64 + - 34 

°' m -0.51 


1290 


T 


Abell 85 


904 


I 


10.460 


-9.303 


0.059 


3.4 


c £C +0.14 

6.65 Zq.u 


1290 


- 


Abell 3571 


4203 


S (VF) 


206.869 


-32.864 


0.039 


3.7 


6 71 + 015 

-0.42 


1300 


F,T 


Abell 3667 


889 


I 


303.129 


-56.841 


0.056 


4.7 


7 30 +0-27 
l3y -0.27 


1380 


T 


Abell 478 


1669 


s 


63.356 


10.466 


0.088 


15.1 


707 +0.18 

-0.18 


1450 




Abell 401 


518 


I 


44.739 


13.578 


0.074 


10.5 


,c +0.43 
°- 3J -0.61 


1490 


T 


Abell 2029 


4977 


s 


227.734 


5.745 


0.077 


3.0 


n 11 +0.20 
y li -0.20 


1570 




Abell 2142 


1196 


s 


239.585 


27.232 


0.091 


4.2 


n 45 +0.65 
y ^ J -0.65 


1610 




Abell 3266 


899 


I(VF) 


67.815 


-61.456 


0.055 


1.6 


9 86 +0 M 
y o " -0.54 


1650 


F,T 



Table 1. Some basic properties of the sample, listed in order of increasing temperature. Coordinates are the centroids used in the X-ray analysis. "Chandra 
observation identifier; ^Denotes either ACIS-I or ACIS-S (observations telemetered in very faint mode are marked with VF); c Galactic HI column, interpolated 
to the cluster centroid using the data of iDickev & LockmarJ <199fj) ; d The mean temperature measured between 0.1-0.2 rsoo (see Section l4~Ti : e An "F" 
indicates that the galactic absorption was frozen at the HI value for the analysis, while "T" marks those clusters for which the projected T(r) was used for the 
3d T(r). All errors are 68% confidence. 



blank-sky data sets . To allow for small variations in the particle 
background level between the blank-sky fields and the target obser- 
vation, we rescaled the effective exposure of the background data 
sets according to the ratio of count rates in the particle-dominated 
7-12 keV energy band. This ratio was calculated for those CCD 
chips not considered part of the main detector (i.e. excluding chips 
0-3 for ACIS-I and chip 7 for ACIS-S observations). Generally 
there is good agreement between the ratios found for different chips 
in a given observation. To avoid the bias caused by contaminating 
point sources in the determination of the background rescale fac- 
tors, we identified an d excluded such features using the it erative 
method described in S anderson. Finoguenov & MohJ 120051) . 

To gauge the sensitivity of our results to variations in the nor- 
malization of the blank sky datasets, we investigated its impact on 
our analysis of Abell 4038 (mean temperature ~3 keV) — the clus- 
ter most susceptible to this effect (i.e. where the background is most 
dominant). We refitted spectra from our annular profile (see Sec- 
tion l4.3> . separately adjusting the background normalization by 10 
per cent higher and lower. In the outermost annulus, this biased the 
recovered temperature by roughly lcr, in the direction of increas- 
ing temperature with higher background. However, the bias rapidly 
diminishes for the inner annuli, dropping below 0.3<r for the 3rd 
innermost spectrum. It is clear, therefore, that our results are not 
sensitive to uncertainties associated with the use of blank sky back- 
ground datasets. 



4 DATA ANALYSIS 

4.1 Cluster Mean Temperature and Fiducial Radius 

The primary focus of this study is on cluster core properties, and we 
have therefore chosen to devote our analysis to the main CCD chips 
for each of the two detectors, i.e. chips 0-3 for ACIS-I and chip 7 
for ACIS-S. For the spectral fitting, weighted response matrix files 
(RMFs) were generated using the CIAO task 'mkacisrmf'. How- 
ever, for some observations this was not possible (as indicated in 
Section[5J, owing to a lack of calibration data suitable for mkacis- 
rmf, and so the older task 'mkrmf ' was used instead. In each case 
the correct gain file was used, appropriate as described in the CIAO 
documentation. 

In a scaling study such as this, it is important to normalize ob- 
servable quantities appropriately, to provide a fair comparison be- 
tween clusters of different sizes. The key parameters of interest are 
mean gas temperature and characteristic radius, and we outline be- 
low a simple scheme in which both these quantities are determined 
self-consistently. Our fiducial radius corresponds to an overdensity 
of 500 with respect to the critical density of the Universe. This ra- 
dius, rsoo, is commonly used for scaling studies and corresponds 
to roughly 2/3 of r2oo I Sanderson & Ponman 2003). 

We defer a full mass profile analysis and direct determina- 
tion of overdensity radii to a later paper, and adopt in this work a 
simple, empirically calibrated proxy for the total mass within rsoo, 
b ased on the mean tempera ture, derived from the M — Tx relation 
of Finoguen ov et alJ 12001). This has the advantage of permitting a 
direct comparison with other observations where poorer data qual- 
ity prevents a full mass analysis being performed. For a cluster of 
redshift, z, the radius is given by l Will is etall2005h 
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Figure 2. A comparison of the mean temperatures from this work and those 
o fllkebe et alj 120021) . showing the line of equality. 



where, 



E(z) = (l + z)Jl + (z + 



A 



(2) 



In calibrating the temperature-radius relation using observed 
clusters, we avoid the bias inherent in equivalent relations derived 
from numerical simulations, which tend to overesti mate the virial 
radii of the coolest haloes ISanderson et al]l2003l) . In using the 
above expression for rsoo we have not allowed for the effects of 
evolution in the overdensity factor (i.e. the '500') with redshift. 
However, for such a low redshift sample, the impact of this evo- 
lution is negligible. 

The mean cluster temperature, T, is measured from the spec- 
trum extracted within 0.1-0.2 rsoo, thus excluding the central re- 
gion where the effects of strong gas cooling can contaminate the 
X-ray emission. The combination of T and rsoo are determined it- 
eratively, by extracting an initial spectrum, finding the best-fit tem- 
perature, then using this to estimate rsoo, and repeating the process 
until convergence is achieved. 

Fig- ID shows a comparison of the mean temperatures derived 
with those from the original ASCA analysis of llkebe et alJ 120021) . 
who used a 2 temperature fit to correct for the effects of any cool 
core. There is good agreement below ~5 keV, but above this our 
temperatures are systematically hotter in almost every case, for 
both CC and non-CC clusters. However, this behaviour is most 
likely the result of differences in the regions from which spec- 
tra were extracted in the two analyses. The most massive clusters 
have large angular sizes, which are significantly larger than the 
Chandra, but not ASCA, field-of-view. Since t he temp erature al- 
most always declines beyond any cool core (e.g. iMarkevitchl 1 998t 
Vikhlinin e tai]|2005t Ipiffaretti et alj f2005). this would lead to a 
higher mean temperature measured with Chandra. 



4.2 Cool Core Clusters 

Having determined the cluster mean temperature, unbiased by the 
effects of central gas cooling, the spectrum of the core region (<0.1 
rsoo) was extracted separately and fitted as before to yield a core 
temperature (see Table[2}- A comparison between this core temper- 
ature and the mean (core-excluded) temperature can then be used to 
quantify the influence of central cooling. By considering the ratio 
of the mean to the core temperature as a discriminator, we define 
cool core (CC) clusters as those systems for which this ratio ex- 
ceeds unity at greater than 3<r significance. This provides a clean 
separation of the sample into 9 CC and 1 1 non-CC clusters, with 
a mean ratio and standard deviation of 1.35/0.15 and 1.0/0.11, re- 
spectively. 

For comparison, in a recent Chandra study Bauer et al. 1 2005) 
find that at least 55 per cent of their sample of 38 X-ray luminous 
clusters show signs of mild cooling, with 34 per cent displaying 
evidence of strong cooling. According to the definition used by 
iBauer et"all all the clusters in our sample could be classified as 
having a cool core (cooling time, t coo \< 10 Gyr). However, their 
sample is more distant (z ~ 0.15-0.4), and they are correspond- 
ingly less able to resolve the inner parts of the ICM where the gas 
cooling time is lowest an d possibly f alls below 10 Gyr in many 
cases. The older study of IPeres et al] 1 19981) reported a cool-core 
fraction of 70 per cent, based on lower-resolution ROSAT observa- 
tions of a complete sample of the 55 brightest clusters in the sky 
in the 2-10 keV band. However, their definition of a cool core is 
different again, requiring that the upper limit to the central t coo i be 
less than the assumed age for the cl uster. 

In any case, both B auer et al] J20051) and IPeres et~aT] fl998) 
base their definition of a cool-core on the gas cooling time which, as 
they demonstrate, is clearly capable of reaching low values (~few 
Gyr) in most — perhaps all — clusters (see also Section lo"2l . By 
casting our definition in terms of a significant temperature decrease 
in the inner 0.1 rsoo, we are able to identify clusters where radia- 
tive cooling has demonstrably impacted the gas properties in a sub- 
stantial way, beyond merely forming at least some gas with a low 
cooling time. 



4.3 Spectral Profiles and Deprojection Analysis 

In order to study the spatial variation of gas temperature, a pro- 
jected temperature profile was obtained for each cluster. Spectra 
were extracted in a series of concentric annuli, centred on the peak 
of the X-ray emission. The radial bins were chosen to enclose a 
fixed number of net cluster counts between 1000 and 3000, depend- 
ing on the quality of the observation. Each spectrum was fitted with 
an absorbed APEC model, as above, to yield the best-fit tempera- 
ture. A characteristic radius was assig ned to each annulus using the 
emission-weighted approximation of McLaughlin ( 1999), 



0.5 



/>/2 
I ' out 



+ r 



3/2 



2/3 



(3) 



To derive estimates of the 3 dimensional gas temperature, we 
used the XSPEC PROICT model to deproject the spectral profiles 
under the assumption of spherical geometry. The deprojection is 
quite slow and susceptible to strong noise fluctuations, so we used 
a smaller number of coarser annular bins — between 10 and 20, 
according to data quality. To stabilize the fitting, the absorbing col- 
umn and gas metallicity we fixed at values obtained by fitting each 
annulus separately prior to the deprojection. For some clusters, it 
was necessary to freeze the absorbing column at the galactic HI 
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Name Mean kT Core kT a kT Ratio 5 flioo CC significance 13 Bautz-Morgan Type Notes' 





(keV) 


(keV) 

Vive V ) 






(kncl 








NGC 5044 


i in +0.03 
115 -0.02 


0.88l°;°° 


1 


ncr + 0.03 

•° -0.03 


4371^ 


11.7* 


- 


C,F 


Abell 262 


9 %A +0.06 
-0.06 


1 OO+0.02 
L - SZ -0.02 


1 


90+0. 04 
• z °-0.04 


6671" 


7.9* 


Ill 


C 


Abell 4038 


J - uz -0.17 


o 9fi +0.14 
°- zu -0.13 





Qn+0.07 
• 9,;i -0.07 


785t! 


-1.1 


III 


- 


Abell 1060 


t ,q +0.07 
• > - ly -0.07 


o 4Q+0.04 





04+O.O2 
a -0.02 




-2.6 


III 


- 


Abell 1367 


1 OQ +0.11 

i - r& -0.10 


q 09 + 0. 10 
°- oz -0.09 





qq+0. 04 
a -0.04 


826^ 


-0.3 


n-m 


M, R 


2A0335+096 


-> A1 +0.14 
3 "> -0.12 


9 q.j+0.03 

z -° -0.03 


1 


40+0. 06 
•*°-0.06 


857+jg 


8.2* 


- 


C, F,M 


Abell 2147 


4 n +0-14 
*- LI -0.14 


4 on+0.37 
^• oz -0.19 





07+0. 06 
°'-0.06 




-2.3 


in 


M 


Abell 2199 


4 36 + - 13 
# - JD -0.11 


qo+0.04 
°- ao -0.04 


1 


-| 1 +0.03 
J± -0.03 


9891^ 


3.4* 


1 


C 


Abell 496 


< 91 +0.58 
Xyl -0.56 


q KC + 0.11 


1 


cc+0.17 
DO -0.17 


II97I70 


3.9* 


1 


F 


Abell 2256 


6.08 +°;f° 


5.64ljj:jg 


1 




1220l3g 


0.7 


ii-iii 


F, M, R 


Abell 1795 


6.30 l^ 7 . 


4.8_l°;«j 


1 


9O+0.04 
• za -0.04 


12461^ 


7.5* 


1 


C, F 


Abell 3558 


6-64 ±° fi 


6.Q4.1°J? 


1 


n +0.09 
lu -0.09 


1288lg 


1.1 


1 


- 


Abell 85 


6-6S ±g;il 


4.80l°;j» 


1 


qo + 0.03 

oo -0.03 


12901" 


11.4* 


1 


F, M, R 


Abell 3571 


-0.42 


_ _„_i_n 17 
'•»O_0.17 





O r + O.04 

S5 -0^04 


f-.^,— -1-1 Q 

1297_^ 


-3.6 


1 


M 


Abell 3667 


7 ™ +0.27 
-0.27 


6 76+ ' 18 
D - 'D-0.42 


1 


qo+0.06 
uy -0.06 


I378I32 


1.5 


in 


F, M, R 


Abell 478 


7 n 7 +0.18 
-0.18 


6 00+ 07 
D - uu -0.07 


1 


qo+0.03 

— 0.03 


14481S 


9.6* 




C,F 


Abell 401 


o +0-43 
°- JJ -0.61 


91 +0.50 

s - zi -o.63 
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Table 2. Detailed properties of the sample, listed in order of increasing temperature. "The core temperature measured within 0.1 rsoo; 'Ratio of mean to 
core temperature; c Number of sigma difference from unity (Cool-core clusters are marked with "*"). ^Clusters with detectable cavities in the ICM are marked 
with "C", cold front clusters are marked with "F", probable merger clusters are marked with "M", clusters with radio relics are marked with "R". All errors 
are 68% confidence. 



value (as detailed in Table0, since unfeasibly low values were oth- 
erwise obtained in many of the annular bins. However, we verified 
that in those bins where the absorbing column was able to be fit- 
ted, the optimum values were fully consistent with the HI inferred 
measurement. Moreover, we found no indication of any radial trend 
in absorbing column in these cases, which might point to problems 
with the calibration. Similarly, in a number of cases the deprojected 
temperature had to be fixed at its projected value, to produce a sta- 
ble fit (denoted by a 'T in the rightmost column of Table0. How- 
ever, these are all non-CC clusters, with approximately isothermal 
temperature profiles, where the smoothing effects of projection are 
minimal, thus any bias introduced by this approximation is likely 
to be small. 



5 NOTES ON INDIVIDUAL CLUSTERS 

In this section, we provide further information about each of the 
clusters in the sample, highlighting key aspects of the analysis spe- 
cific to different datasets. 

5.1 NGC 5044 

This is the only galaxy group in the sample, since a selection based 
on flux biases towards more massive clusters. It is a well-studied 
object, which hosts a cold front and cavities seen in X-ray emission 
iBuote et alJ2003l) . 

5.2 Abell 262 

A known cavity cluster iBlanton et all2004l) . The optical properties 
of this cluster are unusual, and its luminosity function and colour- 
magnitude relation shown signs of contamination from a large num- 



ber of lower mass galaxies, possibly associated with a nearby su- 
percluster (W. Barkhouse, private communication). 

5.3 Abell 4038 

It was necessary to use the projected temperature profile in place of 
the deprojected one for this cluster, in order to stabilize the fitting 
and thereby avoid strong fluctuation between adjacent bins. 

5.4 Abell 1060 

Optical observations in dicate that A 1060 has a dynamically per- 
turbed condensed core iGirardi et al. 1997). However, no obvious 
merger activity is evident from the X-ray emission, so this cluster 
has not been classified as a merger candidate in TableQ The pro- 
jected T(r) data were used for the deprojection, to avoid strong 
fluctuations in the recovered profile. 

5.5 Abell 1367 

This is a complex structured cluster, which is clearly undergoing 
multiple me rging, as seen in both optical (e.g. lCortese et al. 2004) 
and X-rays (Su n & Murray! l_00-l) . The projected T(r) data were 
used for the deprojection, to avoid strong fluctuations in the recov- 
ered profile. 

5.6 2A0335+096 

This is known to b e a cold front cluster 

( Mazzotta, Edae & Markevitch 2003), which also possesses a 
radio-lobe cavity l Birzan et al. 2004). A recent deep XMM-Newton 
analysis of this system indicates th at it is also likely to be 
undergoing a merger with a subcluster l Wer ner etalj2006l) . 
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5.7 Abell 2147 

This cluster is composed of several clumps, with evidence of lumi- 
nosity segregation in the galaxies lLuggeiiri989l) : it has thus been 
labelled as a merger candidate in Table Q It is interesting to note 
that the brightest cluster galaxy for A2147 is not located at the cen- 
tre of the cluster potential (and X-ray peak) l Lugger 1989). 

It was necessary to fix the absorbing column at the galactic HI 
value, since unfeasibly low values were obtained when it was left 
free to vary. The projected T(r) data were used for the deprojec- 
tion, to avoid strong fluctuations in the recovered profile. 



5.8 Abell 2199 

This is a known cavity cluster Ijohnstone et al.l2 002). 



5.9 Abell 496 

This cluster hosts a prominent cold front iDupk e&Whitel2003l) . 



5.10 Abell 2256 

This is a probable c old front cluster in the early stages of merg- 
ing ISun et ai]|2002t) . Our annular spectral bins were centred on 
the main (East) cluster peak, referred to as 'PI' in the analysis of 
ISun et alJj2002T) . 

It was necessary to use the projected T(r) for the deprojected 
T(r), to avoid strong fluctuations in the deprojected profile, aris- 
ing from a complicated 'S' shaped profile. The absorbing column 
was also fixed at the galactic HI value for all bins, to stabilize the 
fitting; there is no evidence from the global spectrum of any signif- 
icant difference between the fitted column and the HI value. It was 
necessary to use the older CIAO task 'mkrmf ' to generate spectral 
responses for this dataset. 



5.11 Abell 1795 

This is a cold front cluster t Markevitch. Vikhlinin & Mazzottal 
1200 ll) . 



5.12 Abell 3558 

It was necessary to use the projected T(r) for the deprojected T(r) 
for this cluster. 



5.13 Abell 85 

A well-known subclump clust er, which also has a col d front and 
is a probable merger candidate iKempner. Sarazin & RickeJ2002l) . 
The prominent Southern subclump and also the Western subclump 
were both excluded from the analysis. 



5.14 Abell 3571 

This c luster is probably in the late stages of merging IVenturi et alJ 
2002) and is located in the Shapley supercluster. It was necessary 
to use the projected T(r) for the deprojected T(r). 



5.15 Abell 3667 

This is both a cold f r ont and merging cluster 
IVikhlinin. Markevitch & Murray! 1200 ll) . It was necessary to 
use the projected T(r) for the deprojected T{r), 

5.16 Abell 478 

This well known cool core clust er hosts a cold front 
iMarkevitch. Vikhlinin & Formanll2003l) a nd possesses pr ominent 
X-ray cavities with coincident radio lobes ISun et alj200l) . 

5.17 Abell 401 

This is a probable merger remnant l Sakel liou & PonmariEooj) . 
The older CIAO task 'mkrmf was used to created a spectral re- 
sponse matrix, owing to lack of calibration data for the newer task 
'mkacisrmf'. It was necessary to use the projected T(r) for the de- 
projected T(r). 

5.18 Abell 2029 

Another well studied cool core cluster which possesses a cold front 
jMarkevitch et all2003l) . 

5.19 Abell 2142 

This is the archetypal 'cold-front' cluster, which likely consists of 
the largely intact core of a poor cluster enclosed within a halo of 
much hotter gas associated with a larger cluster (Markevitch et al. 
2000). This fact accounts for its unusual temperature profile, which 
is atypical of non-CC clusters (see section l6"TY 

The two observations of Abell 2142 (obsid 1 196 & 1296) were 
amongst the very first made by Chandra, and were taken when the 
CCD temperature was only -100C. Since the standard CIAO cali- 
bration does not cover this CCD temperature, these datasets were 
analysed according to the calibration for a CCD temperature of 
— HOC. In addition, the corresponding Markevitch period A blank 
sky background datasets are truncated at 10 keV, so the background 
renormalization was performed in the range 8-10 keV. The analysis 
presented here is based on the 1 196 dataset, but to verify our results 
we also analysed the 1296 observation, and found very good agree- 
ment between the density and temperature profiles. As for A401, 
it was necessary to use the older CIAO task 'mkrmf to generate 
spectral responses for this dataset. 

5.20 Abell 3266 

This is a well-known merging cluster iHenriksen & Tittlevfc o02). 
It was necessary to use the projected T(r) for the deprojected T(r) . 



6 RESULTS 

6.1 Temperature Profiles 

The temperature profiles from the projected annular analysis are 
plotted separately for each cluster in Fig. [3] with vertical dot- 
ted lines indicating the region 0.1-0.2 rsoo used to calculate the 
mean temperature. To clarify the underlying trend in each pro- 
file, the raw data points have been fitted with a locally weighted 
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regress ion in log-log space, u sing a similar method to that out- 
lined in Sanders on et aljl2005l) . This technique effectively smooths 
the data using a quadratic function which is moved along the 
set of points to build up a curve, in an analogous fashion to 
the way a moving average is computed for a time series. The 
algorithm used is implemented in the 'loess' function in ver- 
sion 2.3 of the R PROJECT sta tistical environment package 4 
lR Development Core Tearr]|2006l) . which also provides an esti- 
mate of the la error on the regression, based on the scatter about 
the curve. To compare profiles across the sample, the projected tem- 
perature are plotted as a function of physical radius (in kpc) for CC 
and non-CC clusters in separate panels in Fig. |4] Once again the 
curves represent a locally-weighted regression, this time using the 
R PROJECT 'lowess' function (as used in lSanderson et alj [:2005). 
which provides a larger degree of smoothing. The curves have been 
assigned an arbitrary line style to uniquely identify each system. 

The difference between CC and non-CC clusters is very clear 
from Fig.|4] Within the cool core, there is a remarkable consistency 
in shape of the temperature profile, which follow a power law with 
a logarithmic slope of roughly 0.3. In most clusters this behaviour 
appears to extend down to a few kiloparsecs — well within the cen- 
tral galaxy. The obvious exception to this is the coolest system, the 
NGC 5044 group, which exhibits a sharp break at ~10 kpc, as the 
logarithmic slope flattens within the central galaxy. As this is the 
only galaxy group in the sample, this result may point to a differ- 
ent scaling of hot gas at small mass scales. However, this system 
does display evidence of multiphase gas and morphological distur- 
bances, in the form of X-ray cav ities from AGN activity as well 
as a cold front teuoteetal. 2003). It is therefore possible that this 
flattening in T(r) is an atypical feature. However, we emphasize 
that the possibility of a flattening of T(r) in the inne rmost core 
of clusters (as recently claimed by Donahue et al. 2006, for exam- 
ple) is not ruled out by our data, since the 'lowess' regression tends 
to mask such behaviour, by down-weighting any deviations of the 
innermost few points from the general trend. This disadvantage is 
compensated for by the clarity with which the local regression cap- 
tures the underlying trend in T(r), to highlight the scaling between 
different clusters. 

Fig.|5|shows the variation in projected temperature with scaled 
radius. The size of the cool core region appears to be reasonably 
well approximated by 0.1 rsoo — the radius of the core exclusion re- 
gion used in calculating the mean cluster temperature (Section l4~TY 
However, any turn-over in the temperature profile at larger radius 
is difficult to detect, given the limited field-of-view of the main 
Chandra detectors. This is compounded by the fact that all but one 
of the CC cluster observations analysed here used the ACIS-S de- 
tector, which covers 1/4 of the area of ACIS-I. A further systematic 
difference is evident at smaller radii, where the profiles in Fig. |4] 
extend further in for CC clusters, as a consequence of the bright 
cool core enabling finer annular binning. 

The individual deprojected scaled temperature profiles are 
plotted in Fig. |6| as a function of scaled radius. Here each pro- 
file has been normalized by the mean temperature of the cluster, to 
provide a direct comparison between systems of different mass. It 
is clear that there is a tight spread between the cool core clusters, 
which hold to an approximately power-law form, with a logarith- 
mic slope of ~0.4 (plotted as a solid line with arbitrary normal- 
ization). Conversely, the non cool-core clusters scatter around the 



4 http://www.r-project.org 
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Figure 5. Projected, finely-binned temperature profiles for each of the clus- 
ters, scaled to rsoo and labelled according to the presence of a cool core. 
Each curve represents a locally-weighted fit to the data points, to suppress 
scatter (see text for details). 



locus of isothermality (dotted line), showing no obvious systematic 
variation with radius. 

Motivated by the uniformity in temperature profiles displayed 
by the CC clusters, we have combined all the scaled data points 
for these systems and subjected them to a locally-weighted re- 
gression, to determine a characteristic average temperature profile. 
Each cluster's deprojected temperature profile was scaled by the 
mean cluster temperature (T) and corresponding rsoo • The regres- 
sion was performed in log-log space, using the R PROJECT task 
'loess', in order to provide an estimate of the error on the regres- 
sion, from the scatter about the line. Each point was weighted by 
its inverse variance (computed using the mean asymmetric error), to 
improve the rejection of outliers. The same procedure was applied 
to the non-CC clusters, to provide a comparison. The correspond- 
ing profiles and la error envelopes for the CC and non-CC clusters 
are plotted in Fig.0 The outermost point has been excluded from 
the mean CC and non-CC profiles plotted, in order to minimize any 
bias that this single point can cause at such large radius. 

The behaviour of the profiles at large radii is less certain, 
as indicated by the widening of the regression confidence region, 
resulting from the diminishing number of points. However, there 
is clearly a preference for a peak at ~0.15 rsoo, followed by a 
decli ne, which is broadly consistent with previous an alyses (e.g. 
IVikhlinin et all2005tlpiffaretti et all2005llZhang et alfeOQrj) . 



6.2 Cooling Time 

A simple estimate of the susceptibility of gas to radiative cooling 
can be obtained by considering the timescale over which the gas 
can continue to lose energy at its current rate. Thus defined, the 
cooling time (in seconds) of a parcel of gas with volume, V (in 
cm 3 ), and luminosity, Lx (in erg s _1 ), is given by 
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Figure 3. Projected temperature profiles and lcr error envelopes for each system, ordered by increasing mean temperature from top left to bottom right. The 
curves represent a locally-weighted fit to the data points (in log-log space). Cool-core clusters are denoted by solid lines and non cool-core clusters are plotted 
with dashed lines. Vertical dotted lines bracket the region 0.1-0.2 rsoo, within which the mean temperature is measured. The mean temperature for each 
cluster is indicated, as is the number of annular bins used to construct the curves. 
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Figure 4. Projected, finely-binned temperature profiles for the sample, split by cool core type and using an arbitrary line style to uniquely identify each 
system: the cluster name entries in the legends are ordered by temperature. Both plots are scaled identically in temperature and radius. Each curve represents 
a locally-weighted fit to the data points (in log-log space), to suppress scatter (see text for details). 




Figure 6. Deprojected temperature profiles for CC and non-CC clusters, 
scaled by T and plotted as a function of scaled radius. Each curve represents 
a locally-weighted fit to the data points, to suppress scatter (see text for 
details). 



Figure 7. Mean deprojected temperature profiles for CC and non-CC clus- 
ters, scaled by T and plotted as a function of scaled radius. Each curve 
represents a locally-weighted fit to the data points, to suppress scatter (see 
text for details). The shaded regions represent the \u confidence region of 
the fit. 



t cool = 1.602 x 1CT 9 x 3kTpVflc 



2fiL x 



(4) 



where kT is the deprojected gas temperature (in keV), p is the elec- 
tron number density (in cm -3 ) and the constants /i c and /i are the 



mean mass per electron (1.167) and the mean molecular weight 
(0.593) of the gas. Profiles of cooling time as a function of scaled 
radius are plotted in Fig. [8] colour coded by the mean cluster tem- 
perature. As before the curves represent a locally-weighted regres- 
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Figure 8. Gas cooling time profiles for each of the clusters, scaled to rsoo, 
labelled according to the presence of a cool core, and coloured according to 
the cluster mean temperature. Each curve represents a locally-weighted fit 
to the data points, to suppress scatter (see text for details). 
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Figure 9. Same as Fig.|8] but plotted as a function of radius in kpc. 



sion to the data points. Note that the cooling time in the outermost 
annular bin is unknown (unlike the temperature), since the assumed 
volume is not well defined, which means that the gas density cannot 
be determined. 

For purely bremsstrahlung X-ray emission, t coo i oc \/kT. 
However, the contribution of line emission becomes significant at 
lower temperatures, which weakens this trend. Therefore this sim- 
ple scaling has not been applied to Fig. [8] and the profiles are thus 



largely ordered by mean cluster temperature. Nonetheless, despite 
this spread in temperature, there is a clear consistency between all 
the profiles, with much greater similarity between CC and non- 
CC clusters than is seen in the gas temperature. Not surprisingly, 
the cool core systems have consistently lower cooling times at all 
radii, but they also exhibit a slightly steeper logarithmic slope than 
those of the non-CC clusters. While it is clear that CC clusters reach 
much lower values in the core, it can be seen that within ~0. lrsoo, 
most non-CC clusters have rather low cooling times, falling below 
a Hubble time in all cases. 

It can be seen in Fig. |S| that the locus corresponding to the 
Hubble time intersects the profiles around 0. lrsoo. This is the same 
radius used to excise any cooling-dominated flux from the mean 
cluster temperature estimates (Section fa.U . and also approximately 
corresponds to the point where the CC cluster temperature profiles 
reach a maximum (Fig.0. It is apparent that even in non-CC clus- 
ters, the cooling time of gas within this radius can fall well below 
a Hubble time. However, the non-CC gas cooling time seems to 
reach a minimum level of around 1 Gyr, whereas all the CC pro- 
files fall well below this value. To some extent this behaviour can 
be explained by the fact that cool core clusters permit a much finer 
scale spectral binning, and so their gas properties can be measured 
to smaller radii, where t coo \ is naturally lower. Nevertheless, non- 
CC clusters have noticeably shallower cooling time profiles at all 
radii. 

An alternative perspective is obtained by studying cooling 
time as a function of physical rather than scaled radius, as plot- 
ted in Fig. [5] Here the scatter between the profiles is actually re- 
duced, particularly beyond ~100 kpc. Moreover, the two largest 
outliers in Fig.|8|(the coolest two systems) lie well within the main 
trend, albeit with NGC 5044 showing a significantly steeper outer 
slope. The implication of this is that the cooling time of gas in viri- 
alized systems (either with or without a cool core) obeys a uni- 
versal form across more than two orders of magnitude in radius, 
wi th relatively small scat ter ( a factor of ~3), as found previously 
bv lVoiet & FabiarJ 120041) and lBauer et"ail 12005). This suggests a 
prominent role for radiative cooling in governing cluster proper- 
ties. It is not clear exactly why this should be so, but it represents 
an interesting challenge to theoretical models of cluster properties 
to account for this behaviour. 

Within ~100 kpc t coo i scales roughly as r 1 ' 3 , as indi- 
cated by the so lid black line, which agrees with the findings of 
IVoigt & FabiarJ J2004I) . based on a sample of 16 c lusters with 
strong cool cores. A similar result was obtained bv iBauer et alJ 
(2005), from a sample of 38 X-ray luminous clusters at moder- 
ate redshift. Furthermore, the apparent steepening of the profiles at 
larger radius is also consistent with both these previous studies. 



7 DISCUSSION 

It is clear that clusters can be subdivided into two distinct cate- 
gories according to the presence o r absence of a cool core (e.g. 
iPeres et aJl998tlBauer et all2005l) . How then can we explain the 
mixture of these types in the cluster population? Judging from our 
cooling time results, gas in the cores of all clusters falls well below 
a Hubble time (down to a least ~1 Gyr). This demonstrates that 
the intracluster medium is certainly capable of developing a cool 
core in the inner regions; that it does not always do so implies the 
influence of some counteracting process(es). 

Given the predominantly isothermal temperature profiles 
seen amongst the non-cool core clusters, it is natural to ex- 



© 2006 RAS, MNRAS 000,[T][T4] 



Temperature & Cooling time profiles 1 1 



amine the role of cluster merging in erasing cool cores, which 
would tend to mix gas and smooth out thermal gradients (e.g . 
Fabian. Nulsen & Canizareslll984l : bdge. Stewart & Fabianlll992t 
Buote & Tsailll99a) . The rightmost column of Table [2] indicates 
which clusters show signs of merging activity, and a clear trend 
is apparent: 8 out of 10 'merger' clusters are non-CC, leaving only 
3 out of 11 non-CC clusters showing no obvious signs of merg- 
ing, namely Abell 1060, Abell 4048 and Abell 35 58. Of these 
A 106 has a dynamically perturbed condensed core iGirardi et alJ 
Il997t) and may be a relatively recent merger cluster in which a cool 
core h as not yet been re-established lYamasaki. Ohashi & Furushol 
2002). The CC clusters with indications of merging activity are 
Abell 85 and 2A0335+096. Of these, A85 clearly contains two sub- 
clusters and it appears that t he cool core has surviv ed a merger with 
at least one of these objects l Ke mpner et alj2 002|). A similar event 
may well be occurring in 2A0335+096 IWerner et al]l200d) . thus 
demonstrating that not all mergers are capable of destroying an ex- 
isting cool core. 

If merging is primarily responsible for erasing cool cores, then 
its impact outside this region must be relatively minor. We note that 
there is no evidence that the non-CC clusters are significantly hot- 
ter on average than the CC ones, (e.g. Fig.0, based on the tem- 
perature measured within 0.1-0.2 rsoo (Section |4.H : a compari- 
son of the temperatures of the CC and non-CC clusters, using the 
Kolmogorov-Smirnov Test (implemented in the R PROJECT func- 
tion 'ks.test'), results in a p- value of 0.9195 in favour of the null 
hypothesis that the two sets of data are drawn from the same con- 
tinuous distribution. It can be seen that the cooling time profiles of 
both types of cluster roughly converg e beyond ^O.lrsoo (Fig[8j. 
The head-on merger simulations of Gomez et al. (2002) indicate 
that the residual signature of a cool core may remain for up to 1-2 
Gyr after a disruptive merger, before the onset of turbulence heats 
the gas significantly. Thus it is possible for a cool core cluster to 
have recently experienced a significant disruption, which could ac- 
count for the presence of cool cores and merger signatures in both 
A85 and 2A0335+096. 

Within the core, the apparent 'floor' value of ~1 Gyr seen 
in the cooling time of non-CC clusters could represent the time 
taken to re-form a cool core. If mergers are primarily responsible 
for erasing cool cores, then this timescale can be linked to the char- 
acteristic timescale between such events. Specifically, the fraction 
of non-CC clusters would correspond to the fraction of the time 
between mergers that is required to re-establish a cool core. There- 
fore, if the ~ 1 Gyr 'floor' value represents the time taken to re-form 
a cool core and roughly 50 per cent of the sample possess a cool 
core, this implies a characteristic ti me of ~2 Gyr between mergers. 
By comparison, the simulations of ICohn & White! 120051) indicate 
that, on average, clusters in their sample of 574 have experienced 
at least 4 large mergers since z ~ 2, corresponding to a lookback 
time of ~10 Gyr. This implies a mean timescale of roughly 2.5 Gyr 
between mergers, which is consistent with the above figure. 

Alternatively, cool cores could be disrupted by AGN-driven 
outflows originating in the central galaxy, whi ch is a pla usible 
mechanism for e xplaining their absence (e.g. iMcNamara et alJ 
[2005J: O' Sullivan ^^J20 05|l and is c apable of mixing up gas (e.g 
lO' Sullivan. Vrtilek & Kempner 2005). However, all of the six clus- 
ters in our sample with obvious evidence of AGN-associated cav- 
ities in the ICM (marked with a 'C in the rightmost column of 
Table |2j possess cool cores, reflecting the well established trend 
for cool core clusters to host central radio sources lBurnslll99Cl) . 
On the other hand, AGN may well play an important role in regu- 
lating unchecked cooling in cluster cores, which is well established 



to be significantly suppressed IPeterson & Fabianl2006l and refer- 
ences therein). If AGN activity contributes towards offsetting cool- 
ing in addition to merging, then this would impact the merger rate 
inferred above. Specifically, any contribution from AGN heating 
would require correspondingly fewer mergers, thereby increasing 
the timescale between such events. In that case, the value of ~2 
Gyr deduced above would represent a lower limit to the time be- 
tween major mergers. 

For purely bremsstrahlung X-ray emission the emissivity 
scales as VkTp^ as , which in turn implies t coo i oc VkT. There- 
fore, the apparent universality of the cooling time as a function of 
physical radius (Fig. |9j indicates that hotter clusters must have a 
higher gas density (approaching p gas oc vfeT). However, a given 
physical radius corresponds to a smaller fraction of rsoo in hot- 
ter clusters, which therefore samples gas of higher density, since 
Pgas generally increases towards the centre. Thus, even if gas at 
a given fraction of rsoo had the same density in clusters of all 
masses, there would be some reduction in the scatter in t coo \ plotted 
against physical radius compared to rsoo - Nevertheless, a trend for 
a lower gas density in groups compared to clusters has previously 
been reported (e.g. S anderson et all2003l) . Furthermore, the gener- 
ally good agreement found between the baryon fraction in m assive 
clusters and the the ratio of Ob/^m (e.g. lAllen. Schmidt & Fabian 
l2002tlSanderson & Ponmar]|2003l) . demonstrates that it is the re- 
duced pg as in less massive systems that is anomalous. 

If clusters and groups form self-similarly at any given epoch 
from material of constant density then radiative cooling, or indeed 
non-gravitational heating — operating more effectively in groups 
— could have acted to deplete the gas in the centres of groups, 
thereby lowering the density of the remaining material. However, 
such a process must have occurred subsequent to collapse and viri- 
alization, and so we would expect to observe evidence of systems 
in a intermediate state, which would register as outliers in Fig. [9] 
The lack of significant outliers suggests that the implied variation 
in gas density may already have substantially been in place prior to 
virialization, which implies «on-self-similar accretion. 

Such a scenario is consistent with models in which preheat- 
ing takes places predominantly in filaments, whose linear geometry 
causes the energy injected to impact the density in preference to the 
temperature (as sugg ested by IPonman. Sanderson & Finogueno\l 
l2003UVoit et all2 003). Thus, the gas in the smaller filaments feed- 
ing small haloes could be 'puffed up' more than that in larger 
cluster filaments, lowering the density of material accreted onto 
groups. However, if this picture is correct, it points to an uncom- 
fortable coincidence, namely that preheating of gas prior to accre- 
tion and shock heating produces virialized material of roughly con- 
stant cooling time at any given physical radius. Such an apparent 
conspiracy might only be resolved if some fundamental connection 
between preheating and cooling in haloes could be established. 

7.1 Mergers, Radio Relics and Cold Fronts 

Many clusters are known to possess large-scale, diffuse ra- 
dio sources which are related to the intracluster medium rather 
than AGN iGiovannini & F eretti 2004|). Such sources emit syn- 
chrotron radiation, and are referred to as radio relics. Of the 
clusters in this sample, four are classified as radio relics, 
name ly Abell 85 iGiovannini & Ferettil20q0l). Abell 1367 iGavazzil 
Il978h . Abell 2256 jRottgering et alll994tlClarke & Ensslinl2006T) 
and Abell 3667 (Rottger ing et alJ Il997l) . The non-CC cluster 
Abell 4038 also has a radio halo similar to a relic, but this is most 
likely the remnant of a radio galaxy now located ~18 kpc to the 
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East of the relic fclee et all200lh . Similarly, Abell 2142 has a small 
radio halo whi ch may originate from a sing le galaxy which was ac- 
tive in the past iGiovannini & Ferettil200 0), The generally accepted 
view is that genuine radio relics (i.e. much larger scale than emis- 
sion from a single radio gala xy) are t he product of shock waves 
arising from merging activity lE nsslin et al.ll 9981. and this picture 
is borne out by the fact that all four of these clusters are merger 
candidates. However, A85 also possesses a cool core, demonstrat- 
ing that it is possible for a cool core to survive a merger energetic 
enough to produce a radio relic. 

A related aspect is the prevalence of cold fronts in the sam- 
ple, which are found in 7 of the CC clusters and 3 of the non- 
CC clusters. Of these 10 cold front clusters, 5 show signs of 
merging. Clearly cold fronts are much more common in cool 
core clusters (7/9) than no n-CC clusters (3 / 11), w hich is consis- 
tent with the findings of iMarkevitch et all l2003h . who reported 
that roughly 70 per cent of the nearby cool core clusters in their 
Chandra sample of 37 contained cold fronts. A similar fraction 
has been found mor e recently from an XMM-Newton analysis 
of 62 clusters <Ghizzardi et al) l2005). In any case, the presence 
of cold fronts is certainly not incompatible with the existence 
of a cool core. For very relaxed clusters, like Abell 2029 and 
Abell 1795, cold fronts may result fr om the infall of small sub- 
haloes I Ascasib ar & Markevitchl2 006). rather than being caused by 
more substantial merging. 



7.2 Properties of the Central Galaxy 

Most of the clusters in this sample possess a large central galaxy 
which is coincident with the X-ray peak. The presence of this 
additional potential well is likely to impact the cluster proper- 
ties in the core, so here we briefly explore the characteristics of 
the cen tral galaxy. The Bautz-Morgan (BM) classification scheme 
iBautz & Morgarll970h was devised to catalogue clusters of galax- 
ies according to the contrast between the central galaxy and the 
other galaxies in the cluster (see Bahcall 1977| for a review). There 
are three main categories, which are defined as follows: 

Type I are clusters containing a centrally located cD galaxy 
Type II are clusters where the brightest members are intermedi- 
ate in appearance between cD galaxies and Virgo- or Coma-type 
normal giant ellipticals. 

Type III are clusters containing no dominant galaxies. 

Although the BM type is a subjective quantity, which is sus- 
ceptible to bias l Leir & van den Bergrll977l) . it none the less pro- 
vides a reasonable measure of the evolutionary state of the clus- 
ter. BM classifications are available for 17 of the 20 clusters in the 
sample (the 3 unclassified systems are Abell 478, NGC 5044 and 
2A0335+096) and are listed in Tabled A total of 5 out of the 6 
CC clusters with a BM classification contain type I central galax- 
ies, compared to only 3 out of 1 1 non-CC clusters. Moreover, the 
only CC cluster with a non-type I classification is an anomalous 
system: Abell 262 (type III) exhibits a very unusual colour magni- 
tude relation and luminosity function, with an apparent large excess 
of fainter galaxies, consistent with significant contamination from a 
background supercluster (W. Barkhouse, Private Communication). 

It is clear that BM type I clusters are more evolved than type 
II and III, in the sense that they are unlikely to have experienced 
significant recent disruption. The tendency for CC clusters to be 
BM type I and non-CC clusters to be types II and III is therefore 
consistent with the hypothesis that merger activity removes or oth- 
erwise prevents the formation of a cool core. This behaviour is con- 



sistent with the findings of lBuote & Tsail ll99tjh . who discovered a 
tendency for larger BM types to show evidence of global morpho- 
logical disturbance in the form of a larger power ratio. Similarly, 
iLedlow et alj 12003) find that earlier BM types tend to have higher 
X-ray luminosities, consistent with their hosting cool cores, con- 
firming the tren d for cool core cluster s to have an early BM type 
discovered by i Edg e & StewarJl991 ah . 



While the merger hypothes is has rece ived support from observa- 
tional stud ies (e.g. lAllen et alJl2 00lh as well as earlier s imula- 
tions (e.g. Rick er & SarazirtbOOlk iRitchie & Thomasll2002h . the 
viability of cluster mergin g destroying cool cores has recently 
been called into question ( Gomez et al J 12002: Motl^ydJ 120041; 
IPoole et all 2006). The detailed investigation of lPoole et alj 1 2006) 
indicate that simulations of even head-on, equal-mass mergers are 
incapable of disrupting an existing compact cool core, either by 
heating or mixing of gas. IPoole et al] attribute this new develop- 
ment in part to a better treatment of gas c ooling in t heir simulations, 
which had initially been neglected (Roettiger. Loken & Burns 
ll997tlRicke~ & Sarazin 2001). Furthermore, they incorporate em- 
pirically m otivated initial conditions for their merging clusters 
(outlined in iMcCarthv et"afll2004h . which result in more realis- 
tic and compact progenit or cool cores than those in the study of 
IRitchie & Thomas] 120021) . for example, who nevertheless modelled 
the effects of radiative cooling. 

On the face of it, these new results of IPoole et alj appear to 
rule out the possibility of cool cores ever being disrupted. How- 
ever, this raises the question of why cool cores are not ubiquitous, 
given that the gas in all the non-CC clusters in this sample reaches 
cooling times as low as a few Gyr, although here are a number of 
p ossible explanations for th is (see, for example the recent review 
of Peterson & Fabian 2006, and references therein). Moreover, the 
tendency for non-CC clusters to show signs of recent disruption in 
the form of radio relics and a type 11/111 Bautz-Morgan classifica- 
tion is not otherwise easily understood. Finally, it must be remem- 
bered that radiative cooling and feedback are notoriously difficult 
processes to model in numerical simulations, and that even the most 
sophisticated schemes may fail to adequately capture subtle, small- 
scale physics of the sort that could significantly affect gas dynamics 
in cluster merging. 

Irrespective of any causal relationship between merging activ- 
ity and the presence of a cool core, the distinctive properties of 
the latter reinforce the notion of an intrinsic bi-modality within 
the cluster population. Such a conclusion has important implica- 
tions for the use of scaling relations as a tool for estimating cluster 
properties, such as mass, since any any bi-modality would likely 
dominate the scatter in these relations (see lO'Hara et all2006l for 
example). The extent to which this bi-modality extends to the fun- 
damental property of mass will be addressed in a follow-up paper. 



8 CONCLUSIONS 

We have studied a statistically-selected sample of 20 clusters 
and groups of galaxies, drawn from the flux-limited catalogue of 
Ikebe et al. (2002), using data from the Chandra X-ray satellite. 
The data have been analysed using a non-parametric deprojection 
technique, to estimate the 3 dimensional temperature and density 
structure of the intracluster medium (ICM). We define a quanti- 
tative method to determine the extent to which cooling influences 
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cluster properties: a cool core (CC) cluster is defined as one where 
the ratio of the mean temperature within 0.1-0.2 rsoo to that within 
O.lrsoo exceeds unity at > 3<j significance. Accordingly we find 
that the sample contains 9 CC and 1 1 non-CC clusters. 

We find that there is a clear difference between CC and non- 
CC clusters, with the latter exhibiting somewhat heterogeneous 
properties, although tending to be roughly isothermal within their 
inner regions. CC clusters, by contrast, display a remarkable uni- 
formity in the shape of their inner temperature profiles. The cool- 
ing time profiles display greater uniformity across the sample, with 
non-CC clusters tending to possess longer cooling times at all radii, 
with a slightly flatter logarithmic slope. Nevertheless, even non-CC 
clusters have cooling times much lower than a Hubble time in all 
cases. This fact, together with the high incidence of merger activity 
found amongst the non-CC clusters, indicates that the gas mixing 
and shock heating that this entails may be responsible for erasing 
cool cores or inhibiting their formation. When the gas cooling time 
is plotted as a function of radius in physical units, there is a sur- 
prising decrease in the scatter between different clusters, indicative 
of a universal cooling time profile for gas in collapsed haloes. This 
result suggests that radiative cooling plays a significant role in es- 
tablishing cluster properties. 
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